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Abstract 

The differential cross-sections for processes of tlie type e~^e~ {W~^W~) — >■ luqq are 
determined witli account of background contributions and of anomalous triple gauge 
boson couplings. Analytic expressions for dcr/dsids2dcos6', where 9 is the production 
angle of the W boson, are numerically integrated with the Fortran package GENTLE. 
QED corrections are taken into account in the leading logarithmic approximation. 
The importance of the various contributions is studied for center-of-mass energies of 
190 GeV, 500 GeV, and 1 TeV. 



1 Introduction 



Since the establishment of the electroweak standard model [l]-[3] many precision tests 
confirmed its validity in various respects. One of the poorly investigated features is the 
non-Abelian nature of gauge couplings. W pair production, 

e+e~ W+W-, (1.1) 

provides an excellent way to investigate the triple gauge boson self couplings. 

First calculations of the cross-section for on-shell W pair production in a renormal- 
izable theory, the standard model, were done in the seventies [4, 5]. Already before 
the formulation of the standard model the W width was estimated to yield sizeable 
effects if the W boson is much heavier than the proton [6]. Due to the finite width, W 
bosons decay immediately and the production of four fermions is observed. Production 
of off-shell W pairs, 

e+e- W+W- hfMi, (1-2) 

through the three diagrams of figure 1 was calculated first in [7]. Feynman diagrams 
without an intermediate W pair will also contribute to the four fermion final states 
in (1.2): 

e+e- hf[M- (1.3) 

They constitute the so-called irreducible background and are experimentally not dis- 
tinguishable from the signal diagrams. Classifications and first studies may be found 
in [8, 9], and an overview in [10]. Further, photonic, electroweak, and QCD radia- 
tive corrections must be regarded in order to achieve sufficient accuracy of numerical 
predictions. A huge literature exists on this subject. See e.g. [11]-[21], and references 
therein. 

The properties of triple boson vertices are investigated with different approaches. 
Polarization amplitudes for the most general form of the '-)W^W~ and ZW^W^ ver- 
tices compatible with Lorentz invariance were determined in [22, 23]. Since then, many 
studies appeared on W pair production with anomalous couphngs, see e.g. [24]-[27]. 
For a recent overview, see e.g. [28] . 

The study of the physics of W bosons is one of the main goals of LEP 2 and 
a future high-energy hnear coUider. LEP 2 operates above the W pair production 
threshold at about 161 GeV. Several thousands of W pairs will be produced and precise 
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Figure 1: The doubly resonating CC03 contributions to off-shell W pair production. 
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measurements of mass, width, and couplings of the W boson will become possible. 
Later, at a future linear collider with an energy of 500 GeV or more at high luminosity 
the number of produced W pairs will be even larger than at LEP 2. For a review see 
[29]. 

As mentioned, the process (1.3) may be classified by the final state fermions. In 
this article, we will treat the CCll class, defined by two requirements on the final state 
fermions: (i) they have to belong to two different weak isospin doublets and (ii) no 
electrons nor electron neutrinos are produced. Besides the doubly resonating diagrams 
of figure 1 there are up to eight background diagrams of the types shown in figure 2. 
This depends on the number of neutrinos in the final state: liUil2J^2, h^i<i<i' ■, Q'i^i^2Q'2) 
{li 7^ e). The semi-leptonic CCIO process is of special interest for the study of anomalous 
couplings since its final states offer the most complete kinematical information for an 
experimental analysis of W pair production. 

Present experimental limits on anomalous couplings are not too stringent. In the 
parameter space of aw^i cnw, a-nd the combined limits of LEP and DO are [30] (see 
also [31, 32]) : 

Oiw4> — ~0.03io;o6' 

aw = -0.03^°:°^, (1.4) 
(Xb4> — ~0.051q;2o- 



Here, the identities 



OlWct> = CwSw^Z, (1-5) 
Cw' 



aw = ^7 = — yz, (1-6) 



aB,p ^ Xj - cwSwSz = — — (xz + s^Sz) (1-7) 

are implied, and the anomalous couplings y, are defined in section 3. These 
conventions are in accordance with [28, 33]. 

First signals of anomalous triple gauge boson couplings will be small if any. Since 
the total cross-section is not very sensitive to anomalous couplings, it is advantageous 
to study distributions. 

The semi-analytical expressions of GENTLE for total cross-sections with QED correc- 
tions in the standard model were derived for the signal diagrams in [34, 35] and for the 
background contributions in [36, 37]. With the results presented in this article, GENTLE 
may be used also for predictions of dcr/dcos6', where d is the production angle of one 
of the W bosons. We present analytical expressions for the differential cross-section 
for processes of the CCll class in the standard model in section 2 and the effects of 
anomalous couplings in section 3. In appendices we give some technical details of no- 
tations and the treatment of QED corrections. Numerical results are discussed in the 
corresponding sections. 

The formulae of this article have been implemented in GENTLE version 2 [38] which 
is currently used for experimental studies at LEP 2. 
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Figure 2: Four of the eight singly resonant contributions to off-shell W pair production: 
the d2-diagrams and the U2-diagrams. 



2 The Angular Distribution in the Standard Model 

2.1 The CC03 process 

The CC03 process is defined through reaction (1.2). The fermion pairs fif[ and 72/2 
are the decay products of W~ and W'^: 

W- diu'i, W+ d2u'2, (2.1) 

and have the invariant masses Si and S2- The scattering angle 9 is defined as the angle 
between the electron and the W~ boson. 

The differential cross-section may be written as a sum of s and i-channel contribu- 
tions and of their interference [7] : 

= ^ [dsids2\C'g\s;si,S2,cos9)+C'g'{s;suS2,cos9) 
a cos 9 27rs^ J l 

+C''g'\s;suS2,cos9)]. (2.2) 

We give some notations, including the explicit expressions for the C and Q functions in 
appendix A. 



2.2 Contributions from background diagrams 

We subdivide the background contributions into three parts: 

dcos^ dcos^ dcos^ dcos^ 

The first term contains the interferences between the two s-channel resonant diagrams 
and the eight background diagrams. The second one describes the interferences be- 
tween the t-channel exchange diagram and background, and the third one the pure 
background contributions. 

We denote the various background diagrams by the type of final state fermion 
coupling to the neutral gauge boson. If e.g. an up-type anti-fermion couples to the 
photon or the Z boson, we will call this a -Ui-diagram. In accordance with (2.1), the 
subindex 1 (2) indicates by convention that a fermion of the weak doublet with negative 
(positive) net charge is couphng to the neutral boson. We use the calculational method 
described in [36] and FORM [39]. 
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2.2.1 Background— s-channel interference 



There are 16 interferences between the two ,s-channel signal diagrams and the eight 
background diagrams. Each of these interferences is split up into two products 
and C!"^^!"*: 



da. 



dcos6' 2txs^ 



Jdsids2 Yl E [c;»*^?;«*+cf"^^;!«^] . (2.4) 

i=l,2 a=u,d 



Summation index a stands for up-type or down-type fermions of doublet i. The coef- 
ficient functions are: 



Qsai 



A:,/=7,Z 



(67r2)2 



Rc- 



1 



Dk{s)Dl{s)Dw{si)Dw{s2)D*^{s^_i) 
X Qk [L(e, k)L{e, I) ± R{e, k)R{e, I)] 

X L\Fr, W)L\F2, W)L{f:, l)N,{F^)N,{F2). (2.5) 

The propagators arc defined in (A. 4) and the couphng constants in (A. 5). The two 
independent kinematical functions for the swi-interference are: 

^"v^**^ {2s [s{si + S2) - si- sj £(si; S2, s) (s sif - slj , 

(2.6) 



g'_^'{s,S,,S2) = 



3 1-3 cos2 e 



16 A 

SSS2 



s^siS2 [2ss2C{si; S2, s) s - si S2] 



16 

SSi 



s{si + S2)(l + cos^ 0) + 2siS2 sin^ jC{si] S2, s) 



3S2 



+ - « - 4^2) + ^ s(3Si - S2 - s)(l + cos^ 0) 

o oZ '- 



The logarithm 



2si(si — S2) sin^ 9 
C{s; 81,82) 



A sin^ 9 , , 
H — (si - s- S2). 
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s 



Si - S2 + V\ 



(2.7) 



(2.8) 



Si ~ S2- VX 

arises from integrating the fermion propagators in the background diagrams. 

The ^!"*-functions are proportional to cos^ and, thus, they contribute only to 
the differential cross-section but do not contribute to the total cross-section. After 
integration over cos^^, (2.7) yields eq. (3.1) of [36]. 

One may obtain the SM2-interferencc by exchanging Si and 82 in the stii-interference: 



g±'{s,8i,82)^gr{s,s2,s^). 



(2.9) 



To construct the kinematical functions with the down-type fermion coupling to the 
neutral vector boson, one may use the symmetry: 



Qt'is,Sl,S2) ^QtHs,S2,Si) ^^g'^'i8,8i,82). 



(2.10) 



The coefficients of the V violating contributions in (2.4), vanish for pure photon 
exchange. 
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2.2.2 Background— t-channel interference 

The i-channel background interference is: 



d cos 6 27rs^ 



i=l,2 a=u,d 



(2.11) 



Due to tiie neutrino exchange in the t-channel, only left-lianded particles contribute 
and, therefore, only one combination of couplings appears: 



rRe- 



(67r2)2 Dw{Sl)Dw{s2)Dl{s)D*^{ss-^J 

X L\E, W)L{e, k)L{fl k)L\F,, W)L\F2, W)N,{F,)N,{F,). (2.12) 

The kinematical functions are exceptionally asymmetric since the integration over 
the fermion propagator in the background diagrams is performed, while the neutrino 
propagator (A. 10) in the t-channel diagram is still present. The kinematical function 
for the iiii-interference is: 



(S,Si,S2j = 

— 1 I 3 cos 



A [4 

, sin^ 9 
+ A 



— {s + si- S2) + 2sjC(si; S2, s) 



[2siS2{s2 — Si) - 6s^S2(si + S2)C{si; S2, s) - 3ss2{s + S2)] 



+ 



16 



[{s - sif - si] + 



SSi 



+ 



SSIS2 



-ss2>C(si; S2, s) (5s sin^ 9 + 4si 



+ 4^2) - l{3sl - 2ssi + 4siS2 - 7sl + 30ss2 + Qs^) sin^ 6 - -{3sl - 2sl 
8 2 



-S1S2 + 2ssi) 



3s S2 

H — C{si; S2, s)([4siS2 + sl + si — s{si + S2)] sin^ 9 



A r 9 9 / \n\ Sill ^ / n *? o *? 

- 4[SiS2 + Si + S2 - s[Si + S2)J) H (2SiS2 - OSi + 3S2 

o 

- 14ssi - 3s^) +^(5siS2 - 2siS2 - 3s2 + 5ssiS2 + 3s^S2)| . (2.13) 
The expression for the tdi-interference becomes quite compact using (2.13): 

^?*''Hs,Si,S2) = -^?*"^(s,Si,S2) 



3SS2 



A 



sin^ 9 
4:U 



{(s + Si + S2)[Si(2Si - s - S2) - (s - S2) 



-2 

+s 



SSi(s - Si)^ + SS2(s - S2)^ + SiS2(si - S2)^ >C(si; S2, s) | 



s{si + S2) - si - si C{si; S2, s) + 7: (s + -Si)^ - st 



.(2.14) 



The integral over cos^ of (2.13) yields ^ccii of (2-14) yields Gcci'i^ (eq- (3.12) 
in [36]). 
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The remaining two Q functions are easily constructed: 

G''''(s,s^,S2)^G''''(s,S2,s^) (2.15) 

and 

g"''{s,Si,S2) = g"''{s,S2,Si). (2.16) 

2.2.3 Pure background 

The pure background contribution is: 

^^^^f^^^^^^ ^ ^ Icf^gf^ +C^b^g^b^] (2.17) 

Cl cos f/ ^TTS J 1 J „■ „■ 1 o 

a,o=u,a i,j=L,2 

Again, we have to introduce additional coefficient functions C_ compared to the 
total cross-section, where only C+ functions appear: 

^ khz (6^')' D,is)Dt{s)D^{s,^,)D*^{ss^^) 

X [L(e, k)L{e, I) ± R{e, k)R{e, I)] 

xL{ft,k)L{f^,l). (2.18) 

The potentially 2 x 64 kinematical functions in (2.17) can be reduced to 2 x 16 
functions in a first step since the 7 and Z exchange differ only in the coefficient func- 
tions (2.18). With 

the number of independent ^-functions is further reduced to 2 x 10. 

Finally, we will need only five kinematical functions to express them all. 

The simplest cases are the squares of the various background diagrams {a — b and 
i — j in (2.17)); they are given by: 



/!!'"'(5>5i,S2) = -^^^SS2 |->C(Si;S2,s) {s-Sif-sl +S-Si-S2^ 



(2.20) 

and 

e?rH«,«i,«2) = l^^^^^^ss^sl[C{s^■, S2,s){s2-s, + s) + 2] + cos' e) 

3 

+ ^^■^2(1 + cos^ 9) [sC{si; S2, s){si - S2 - s) - 2s - Si] 

+ ^si{s + 3s2). (2.21) 

By integrating (2.21) over cos^ one gets Qccii (3-3) of [36] while (2.20) vanishes. 
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The other interferences between background diagrams of the same doublet are: 

Gr'(s,s,,S2) = e?rn«,«2,si), (2.22) 



(2.23) 



gi'^^'is^s.^s^) = ±gr'{s,S2,s^). (2.24) 

With the aid of a neutral current function, one may prove the relation: 

ST^'is^s.^s^) = gi"''{s,Si,S2)+gi''''{s,Si,S2)-ss2gg^{s,s,,S2). (2.25) 

The function ^i^f* (s, Si, S2) may be found in appendix C. 

The functions S"^^^ vanish in the neutral current case and the analogue of (2.25) is: 



G'^'''' {s, s^, S2) = - [g'j''' {s, s,, S2) + gi''^' (s, Si, S2] 
The expressions for the other doublet are: 

QI'''\S,S^,S2) = gi"''{s,S2,S^). 



= 0. 



(2.26) 



(2.27) 



The interferences between diagrams from different doublets are more complicated. 
Here, we have 



and the lengthy expressions 



slC{si] S2, s) - slC{s2; s, si) + 



S2 — Si 



2 I 2 
S1 + S2 



(2.28) 



■'Uid2 I 



s, 81,82) = -18 



*^'®l'^2/i I „;„2 /)\ 2 



-3s 



X 



_ ^1 + sin 6)s SiS2C{si;s2,s)C(s2;8,8i) 

Si£(s2; 8, 81) + slC{si; S2, s) 
sin^^ scos^6', , s'^ 8182(1 + siii^ 9) 

+ — r^(« - ^) + 

s(si - 82 



X 



+ 



+ 



4A ' 2A2 

S1>C(S2; 8, 81) - 8lC{si; 82, S) 

3ssiS2(l + sin^ i 



2-3s- 



8 — 3a 



X 



8A + 
3s2siS2(l + sin^6') 

4A2 
s(l + cos^ 6) 
16A 

sin^6' r24ssiS2(2s - a) 



2A2 



l-3s 



8 + a 
A 



222 12ssiS2(s — c) 

8 ~ Si — 82 



A 



32 



4sSiS2 + 3(s? + 8l) - (3s^ + 7siS2)(J 



A 



— 82 — IOS1S2 



(2.29) 
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and 



scos9 J 27s.s{s2 



2A2 



s{a - s)C{si; sa, s)C{s2; s, si) 



+ (si-S- S2)C{si; S2, s) + {s2-s- si)£(s2; s, Si) - 2 
9sSiS2 



+ 
+ 
+ 
+ 



2A 



s[3siS2 + s{a - s)]C{si; S2, s)C{s2; s, Si) 



s{si - s) + S2 [Si - S2 - 



S2KS2 - S- Si) - -SSi 

3si 



£(si;s2,s) 

5 



C{s2]s,si) - + 



6s^S2>C(si; S2, s)C{s2] s, si) + 3ss2JC{si; sg, s) 



-s(3s2 + 2si)>C(s2; s, Si) - (^s + si + ^sa) | ■ 
In (2.29) and (2.30) we use the abbreviation 

(7 = Si + S2- 

Further, 



gi''''{s,Si,S2) = -SSiS2g^,i^{cOSe,S,Si,S2)-gi''''{s,Si,S2) 



(2.30) 



(2.31) 



(2.32) 



The neutral current function ^^f'(cos s, Si, S2) can be found in appendix C. The 
integral of (2.29) is ^ccii defined in (3.10) of [36] and that of (2.28) and (2.30) vanish. 
The remaining kinematical functions are 



and 



^?f"^(s,Sl,S2) = gi"''{s,S2,Si) 
gt'''{s,Si,S2) ^±gi'''%S,Si,S2). 



(2.33) 
(2.34) 



2.3 Numerical results 

Numerical results are obtained with the Fortran program GENTLE [38], version 2.02. 
QED initial state radiation (ISR) is treated as described in appendix B and in [38]. 

Wc use the numerical default input values, e.g. Mw = 80.230 GcV, Tw = 2.0855 
GeV, Mz = 91.1888 GeV, Tz = 2.4974 GeV, sin^ = 0.22591, aem = 1/137.0359895, 
as = 0.12, no Cabibbo mixing, and the GENTLE flag settings 

IPROC = IINPT = lONSHL = IZETTA = 1 

ICONVL = IIQCD = IOCS = IMAP = IRSTP = IMMIN = IMMAX = 1 
IGAMWS = IGAMZS = 1 

IGAMW = ITNONU = IQEDHS = ICOLMB = IZERO = IBIN = IRMAX = 0. 

The flags IBORNF, IBCKGR, ICHNNL are varied in an obvious way. For calculations within 
the standard model lANO is set equal to 0. The flags IGAMWS and IGAMZS are chosen 
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COS0 

Figure 3: The ratios CC09/CC03, CC10/CC03, and CC11/CC03 without QED corrections. 

such that the boson widths are taken to be constant. For the related problems with 
gauge invariance see [40, 41, 10]. 

In figure 3 the net size of the background effects at a center-of-mass energy of 
500 GeV is shown as the ratio of the signal plus background cross-section to the signal 
cross-section 

^^ d.ccn/dcos^ _ (2.35) 
do-cco3/dcos6' 

The background contributions may become sizeable for large scattering angles. For 
extreme backward production, cos^ = —1, the effect is larger than 30%. At = 190 
GeV, dcTft/d cos 6 (see (2.3) is less than 0.3% of ctccos in the whole range of the scattering 
angle; for more details see [42]. 

In tables 1 and 2 we present numerical data which may be of some use for preci- 
sion comparisons of different numerical programs. For this purpose, we ran GENTLE 
at high numerical precision but still a reasonable computing time at a PC with a 
Pentium 133 MHz processor. The numerical reliability was controlled by varying the 
parameters e (for the relative error of Simpson integration) and 6 (a technical cut pa- 
rameter improving numerical stability in some edges of the phase space) in GENTLE. 
The technical uncertainties in the last digits shown in the tables are of the order one 
or smaller. Without ISR the cross-sections are obtained with e = 10~^ and 6 — 10~^. 
With ISR, the corresponding values are e = 3 x 10~^ and S = 10~^ (if IBCKGR=1 and 
IB0RNF=1, then S = lO'^). 
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COS0 

Figure 4: Differential cross-section for e'^e" n'v^ud with various corrections. 

Although there are different couphng constants and even a different number of 
contributing diagrams for the various CC09, CCIO, CCll processes, one observes only 
small differences between the various cross-sections, especially at LEP 2 energies. One 
may suspect that this is due to certain relations between the relevant coupling con- 
stant combinations in the cross-sections which make the latter being dependent on 
only weak iso-spins. For example, the s-channel background interference contributions 
Qsuigsui j^QsdiQsdi identical for all processes of the CCll class. Here, only the parity 
violating terms lead to different background effects in this interference. Similarly, for 
the t-channel background interference the flavour depending combination t/*"* + is 
suppressed by cancellations as it can be seen in (2.14). Obviously, the effects which are 
similar for all final states of the CCll class are the numerically dominating background 
corrections. 

QED corrections are shown in figure 4 at ^/s = 500 GcV. There is a considerable 
cross-section enhancement for cos^ < —0.5. The background effect in this region of 
cos 6* is reduced, but still sizeable. 
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x/i (GeV) 




C'CCOS l,puj 


C'CC09 W^j 


C'CCIO l,puj 


C'CCll l,puj 


190 


U.O 


u.uy^oouo 


u.uouo'yiOO 




n 9QzL'^fin 

u.zy^oou 



u.u 


n 91fi7Q1 


n nfiQ74Q7 


n 917941 


u.u / uuuu 


U.O 


n 7Qn4'?9 




n 7Q1 ni 9 


9 zLfi'^R'^ 

.^UOOtJ 


500 


U.O 






U.UUUOc/ZiUO 


u.u lUOOOU 


0.0 


0.0215908 


0.00654501 


0.0204228 


0.0637232 


0.8 


0.212368 


0.0682179 


0.212503 


0.661957 


1000 


0.8 


0.00211188 


0.000103725 


0.00128268 


0.00107275 


0.0 


0.00614750 


0.00155423 


0.00484706 


0.0151154 


0.8 


0.0535172 


0.0170530 


0.0530947 


0.165310 



Tabic 1: Differential cross-sections without ISR. The CC03 cross-section is calculated 
with the branching ratios for the CCIO process. 



^fs (GeV) 


cos 6 


o^ccos (pb) 


0-0009 (pb) 


o-ccio (pb) 


o-ccii (pb) 


190 


-0.8 


0.09064 


0.02915 


0.09077 


0.2827 


0.0 


0.1971 


0.06340 


0.1975 


0.6150 


0.8 


0.6868 


0.2206 


0.6871 


2.140 


500 


-0.8 


0.01248 


0.003749 


0.01170 


0.0365 


0.0 


0.02590 


0.008036 


0.02506 


0.0782 


0.8 


0.2308 


0.07418 


0.2311 


0.7198 


1000 


-0.8 


0.005282 


0.001441 


0.004510 


0.01412 


0.0 


0.007757 


0.00223 


0.006939 


0.02163 


0.8 


0.0613 


0.01963 


0.06114 


0.1904 



Table 2: Differential cross-sections with ISR. The CC03 cross-section is calculated with 
the branching ratios for the CCIO process. 



3 Anomalous Couplings 

We now extend the Lagrangian of the standard model by anomalous triple gauge bo- 
son couplings. We allow terms that obey Lorentz invariance and CV invariance. In 
addition, for the electromagnetic interaction we forbid C or V violation and will not 
modify its strength. 

These conditions are fulfilled by the Lagrangian proposed in [43] : 



C — —ie 

— ie cot 

— ieSz 

— ie 



'Ml 



Zf, (w-^"'w+ - w+f^w-) + Zi.^w+^w-" 



iex^Ff.^W+^W-' 



- lexzZ.^W^m''' 



+ ^daZpa faw-"w^+" - + d^w+^w-" - d^w+'^w-'') .(3.11 
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The term with the dual field tensor Z violates both C and V. With a multipole 
expansion, one gets the electromagnetic charge Qw, the magnetic dipole moment \iw^ 
and the electric quadrupole moment qw [44, 28]: 



Qw 
qw 



e, 



2Mw 

e 



Mi 



+ + 

2-(l + ^7 -2/7)- 



(3.2) 
(3.3) 

(3.4) 



w 



The anomalous couplings x-y, xz, y^, yz, zz, and 5z produce additional contribu- 
tions to the cross-section of W pair production. The largest contributions will come 
from resonant diagrams (section 3.1), but others are also coming from the interference 
between anomalous s-channel signal diagrams and background (section 3.2). 



3.1 Anomalous contributions to the CC03 process 

We write the cross-section for doubly resonant scattering with anomalous couplings in 
the following form: 



J ano 
'^"CC03 

dcos6' 



27rs2 



dsids2 



Y.^nmQnm{s] Si, S2, 0086) + E<^n ^^n («; ^2, COS 6) 



(3.5) 



The sums over n,m run over x, y, S, and z and the standard model couplings. The first 
sum in (3.5) describes the s-channel interferences and the second sum the anomalous 
■st-interferences. The coefficient functions are: 



nm 



k,hz (6^')' \Dw(siW\Dw{s2WD,{s)Dl{s) 
X 9l9T\'^ + (1 - 

X L2(Fi,iy)L2(F2,PF)7Ve(Fi)7V,(F2), 



E 



(67r2)2 



Re- 



1 



X 5«L(e, V)L\F^, W)L\F2, W)L^{e, W)N,{F,)N,{F2), 



with 



and the standard model couplings 



5?^ 



gSW: 



gcw 



(3.6) 



(3.7) 



9^ 




gz 


= gswxz, 






gswv^ 
Mi ' 


gl 


_ gswvz 


(3.8) 


g'z 


= gsw^z, 


g'z 


_ gswzz 





(3.9) 
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The constant is defined as follows: 



KT-Kf = L{e,k)L{e,l) - R{e,k)R{e,l) for m ^ ^ 
Al"^ = L{e,k)L{e,l)+R{e,k)R{e,l) otherwise 

The 61 in (3.6) is the Kronecker symbol. Note that the pure standard model con- 
tributions are already treated in section 2.1 and should not be counted twice. 
The anomalous kinematic functions in eq. (3.5) are for the st-interference: 



Si' = -,s 

1 



(Sl + S2) (^S - Si - S2 



+ - sin^ e 
U J 4 



-SSIS2 



S — Si — S2 — 



2SiS2 



5'' = h'' 



sin^ 9 

2(si + S2) — (si(s - Si) + S2{s - S2)) 



while for the s-channel contributions: 

1 



xy 



^sx — ^xS 

QS 



128 



:XS 



Si + S2) (1 + cos^ 6) + s sin^ 6 



—XssiS2{l + cos^ 6), 



1 

128' 



As 



4(si + S2) + {s — si — S2) sin^ 9 



-——\ss1S2 2ssin^ ^ + (si + S2)(l + cos^ ^) 
128 '- 



7^AssiS2, 
16 

-^\^s{si + S2){l + cos^ 9), 



»S(5 — ^55 



32 



A' 



2s(si + S2) + I 3siS2 + ^ J sin^ 9 



^xz 



'yz 



1 3 

— A2s(si + S2) cos 61, 
64 

^A5ssiS2Cos6l, 
-^A5s(si + S2)cos6l. 
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3.2 Interferences of anomalous contributions with background 

Finally, we treat interferences of the anomalous s-channel diagrams with background: 

I asias2 

a=u,di=l,2 n 



J^ano 



d COS 9 27rs2 
The coefficient functions are 



/ ds.ds, + c'-:^Q-:^\ . (3.25) 



Csai 



Re- 



1 



(67r2)2 Dk{s)D^{s)Dw{si)Dw{s2)D*y^{s^. 
X gUHe,k)L{e,l) ± R{e,k)R{e,l)] 
X L'{F^, W)L'{F2, W)L{f:, l)N,{F^)N,{F2), 



(3.26) 



where stands for all couplings given in eq. (3.8). 
The anomalous kinematical functions are: 



gTU^,s^,S2) = 



3 cos 6' 
32^ 

3 1-3 cos^ 9 



SS2 |2s[s(si + S2) - sj- sl]L{si; S2, s) + (s + si)^ - sjj , 

(3.27) 

s^siS2 [2ss2>C(si; S2, s) + s - si + S2] 



32 A 

3s 

- — (1 + cos^ 9)ss2{si + S2)C{si; S2, s) 
+ ^{l-3cos^9){s + S2-Si) 



+ 



si - si - s{si + S2) - 



\sm'^9 



(3.28) 



Gt:'Js,s,,S2) 



G'+Zi'^S,,S2) 



16 

3 cos^ f . , 2 

— —-^SSiS2 |2ss2[2s - (Si + S2)\C{Si; S2, s) - 28^ + 2SiS2 

+ 3ss2 - ssi + , (3.29) 
SS1S2 f ^ 1 — 3cos^ 9 



64 



6- 



A 



-SS2 {s[s — (Si + S2)]C{si; S2, s) + S + Si - S2} 



+ (1 — 3cos^ 9)[s — 2ss2j0,{si; S2, s)] — 16ss2jC{si; S2, s) 

- 8(S -Si + S2)> , 



(3.30) 



^!!^^(s,si,S2) = ~^^s\^6ssiS2[2ss2C{si;s2,s) + S - Si + S2] 

+X 6ss2(si + S2)jC(si; S2, s) + s(2si + 3^2) - - Ss^ + 3^2 } , 

(3.31) 

ST,li^,Si,S2) = —{l + cos^9)ss2{2s[sl + sl-s{si + S2)]C{si;s2,s) 

+sl-{s + sif}, (3.32) 
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g'_''}{s,s,,S2) = gi^Ks,si,S2), (3.33) 
g'+^}{s,s,,S2) = g'+^'{s,s,,S2). (3.34) 

The remaining kinematical functions can be calculated with the equations 

g^J^lis, S2) = g£i{s, S2, s,) = ±^-^(5, S2, s,) = ±^-^(5, s,, s^), (3.35) 

gt:As.SuS2) ^gt^{s,S2,S,) ^Tg'^;^,{s,S2,S,) ^Tg'^^l{s,S,,S2), (3.36) 

where a stands for x-y, xz, y-y, yz, and 5z- 
3.3 Numerical results 

In figure 5 we show the bin-integrated differential cross-section for all the six anomalous 
couphngs at 190 GeV. In each case only one anomalous coupling is allowed to differ from 
zero. The figure is in excellent agreement with an analogous figure in [43]. Comparisons 
with the Monte Carlo event generator WOPPER [45, 46] show also agreement within the 
statistical accuracy of the MC program. 

As an application, we shortly describe a study on the discriminative power of W 
pair production with respect to parity conserving and violating anomalous triple boson 
couphngs. At ^/s = 500 GeV with an integrated luminosity £ = 50 fb"-*^, about 80 000 
semi-leptonic W pair decays are produced. The anomalous couplings appear in the 
cross-section at most bilinearly. Allowing e.g. for two anomalous couplings A and B 
simultaneously, one may use the ansatz: 

CTtheor = (7^^ + A(7i + ^Vn + S(72 + + AS(7i2. (3.37) 

After having calculated cr^'^, (Ji, fTn, ...with GENTLE (or another program) within a 
given model and for definite experimental conditions, one may confront experimental 
data with the predictions. For a study of sensitivities, we use crtheor for the simulation 
of <^meas ± \J(^meas/ (S-C), the assumcd measured cross-section with la deviations of the 
counting rates for the sum of all six semi-leptonic production channels. For definiteness 
we use for Umeas the standard model prediction cr^^ and apply no experimental cuts. 
The solutions of eq. (3.37) for A and B are eUipses in the plane. Allowed pairs of 
coupling values are located in the area between the two limiting ellipses. For the 
sample analysis, we use two obscrvables: ap and 0"^, the forward and backward cross- 
sections. The forward (backward) cross-section is defined by the requirement that 
the angle between the momenta of the e~ and the W~ is less (more) than 90°. We 
choose these observables since they may be used to form the total cross-section a tot — 
ap + as (arising from cross-section parts even in the production angle) and the forward 
backward asymmetry Aps = {(t p — a b) I c^tot (arising from odd cross-section parts). For 
two different sets of anomalous couplings, the two rings with allowed values derived 
from (Tp and (Tb overlap almost totally in the case of V conserving couplings x^, 5z- 
When replacing 5z by the V violating coupling zz, the allowed ranges overlap much 
less since the forward-backward asymmetry is more sensitive to this coupling. All this 
is nicely seen in figure 6 (where we also took ISR into account). There, one further 
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Figure 5: The ratio of cross-sections with anomalous couplings to standard model cross- 
sections without background and without ISR corrections at y/s — 190 GeV. In each 
figure only one anomalous coupling differs from zero. 
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GENTLE 



GENTLE 




may notice that as is more sensitive to anomalous couplings than ap although the 
relative statistical error of the latter is much smaller. This is in accordance with the 
properties of the angular distributions in figure 5 (see also figure 2 in [42] for a center- 
of-mass energy of 500 GeV). A similar discussion has been performed for other pairs 
of anomalous couplings in [47]. 



4 Summary 

We determined semi-analytical background and anomalous contributions to the differ- 
ential cross-section for W pair production in processes of the CCll class. With the W 
production angle as an additional parameter, the expressions are not as compact as 
those for the total cross-section. 

By performing numerical calculations with the GENTLE package we illustrated the 
effects of the background contributions. At energies of about 500 GeV or more, back- 
ground has sizeable effects especially for backward scattering. At energies of about 
190 GeV background is less than 0.3% for all scattering angles. Gontributions from 
anomalous couplings arc strongest in the same region. Therefore, they cannot be stud- 
ied without taking background properly into account. 

The present calculation has several limitations, of more importance at higher en- 
ergies: the virtual corrections are not taken into account; for angular distributions, 
the QED ISR radiator function used is only an approximation for the real photonic 
corrections (both with respect to the 0{a) part and to higher order corrections); the 
treatment of finite width effects may be refined. One should estimate the net effect to 
be of the order of up to few per cent at s/s = 500 GeV. Thus, GENTLE is certainly an 
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appropriate tool for the study of W pair production at LEP 2 while at higher energy 
it may serve as a playing ground for sensitivity studies but may not replace a more 
complete calculation. 
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A The CC03 Process 

The coefficient functions used in (2.2) are: 

' .Re. ' 



(67r2)2 \Dwis,W\Dw{s2W 

X L\e, W))L\F,, W)L\F2, W)N,{F^)N,{F2), (A.l) 

,£;;^(67r2)2 \Dw{s,)\^\Dwis2WD,is) 

X gkL{e, l)L\Fi, W)L\F2, W)L\e, W)N,{F,)N,{F2), (A.2) 
2 



k,hz (6^')' \Dw{siW\D^{s2)m{s)Dt{s) 
X gkgi[L{e, k)L{e, I) + i?(e, k)R{e, I)] 

X L2(Fi, W)L\F2, W)N,{F,)N,{F2). (A.3) 
The denominators of the boson propagators are: 

Dvis) ^s-M^ + iMyVv, (A.4) 
and the coupling constants in the standard model are: 

g-r = gsw = e, 9z ^ gcw, 



L{f,W) = R{f,W) = 0, 



(A.5) 



We use Qe — —I and /| = — |. The colour factor Nc is 1 for leptons and 3 for quarks. 
For the kinematical functions Q we quote the expressions from [36] : 



e' = i 



^ , , A . 2 ^ XsiS2 sin^ 9 
2s(si + S2) + ^ sm^ 9 + 



(A.6) 
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7 St 



{s - Si - S2) (2s{si + 52) + ^ sin^ 9 



S1S2 
17 



(4s(si + S2) - Asin^ 6*) 



2s(si + S2) + (3S1S2 + ^ j sin^ 9 



with 



X = s'^ + si + si — 2ssi — 2ss2 — 2siS2 
and the denominator of the neutrino propagator ti,: 



— ^ (^s — Si — S2 — y/Xcos9^ . 



(A.7) 
(A.8) 

(A.9) 
(A.IO) 



B QED Corrections 

The differential cross-sections are calculated in the rest system S' of the W boson pairs 
(or, equivalently, of the final state fermion pairs). If energetic photons are radiated from 
the initial state, S' differs from the laboratory system S where the production angles 
are determined experimentally. The corresponding Lorentz boost will be described 
in appendix B.2. An emission of photons from e~ or e"^ leads to different relations 
between the W production angle in E and in E'. Thus, we have to use the structure 
function approach for a description of ISR since here the energy loss of each initial 
state particle is known. 



B.l Structure function approach 

In the structure function approach [48, 10], the initial state photonic corrections are 
taken into account by convoluting the tree-level cross-section twice with the structure 
function D{x, s) (with the structure functions as described in section 2.3 of [38] and 
references therein): 



do-QED(g) 

dsids2dcos6' 



- ■ 7 = 1 2 

ninin -.mm " ■^'■^ 



d cos 9i 



dcos6' 



da{xiX2S, si, S2) 
d cos 9!i 



with 9'i = 9'i{s, Si,S2,Xi,X2, 9) and the lower integration boundaries 



> 



XiS 



(B.l) 

(B.2) 
(B.3) 
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The sum in (B.l) indicates that no, one, or two solutions may exist for 9[ (defined in 
E') at given values of the parameters in E. The Jacobean is easily derived from (B.23): 



d cos 9[ 



(B.4) 



dcos6' 



PI2 + ^2(1 - (51^ sin^ 9) - 2v(3i,2 cos 9 

r , 1 - 6^ 

X /5i,2 - f cos6^ ± - cos^) — - — - 



-13/2 



V 



(B.5) 



b l±vb cos 9 



The b is given in (B.22), the velocity v of the W'^W~ system in the laboratory frame 
E in (B.12), and the velocities /3i^2 of the W~ in E in (B.17). For doubly resonant 
diagrams the cross-section has to be multiplied by the Coulomb correction C{xiX2s) 
[49, 50, 51]: wc follow [50] as described in [36]. 

For applications and comparisons with Monte Carlo programs [52, 53] it might be 
more convenient to determine not the differential cross-section itself but to perform a 
bin- wise integration: 



Such an integration may be trivially performed analytically in E' in view of the rela- 
tively simple angular dependencies and computer time may be saved. Of course, the 
boosted integration boundaries have to be determined. For a given angular bin in 
the laboratory system, there may exist zero, one, or two bins to be integrated over 
in the boosted frame. More details on this may be found in section 2.4 of [38]. The 
bin-integrated cross-sections are used in GENTLE for CC03 processes. 

Finally, a remark on the use of the structure function D{x, s) might be necessary. 
This structure function is determined for the total cross-section only. Thus, for the 
differential cross-section it has to be considered as an approximation. 

B.2 Lorentz boost 

We will denote 4-momenta in E' as p' and in E as p. In E, the momenta of electron 
and positron are: 



E = v^/2 denotes the beam energy. In E', the sum of the spatial momenta of the two 
particles vanishes and one gets in this frame: 




(B.6) 



Pe+ = EX2{1,0,0,-1). 



(B.7) 
(B.8) 



p;. = 0,0,1), 



(B.9) 
(B.IO) 



p',+ = £;v^(i, 0,0,-1). 



Applying the transformation formula 



P3 - vPo 



(B.ll) 
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on one of the beam particles, one may derive the relative velocity of the two Lorentz 
frames 

Xi - X2 
V = . 

Xi + X2 

In E', one may choose the momenta of the W bosons as follows: 



P'w- 



A 




A' 



Pw+ = \ 7:7 + «2, -VT-ysin^', 0, -W — cos^' 



where s' — AxiX2E'^ is the reduced center-of-mass energy and 

A' = X{s', 81,82). 



(B.12) 

(B.13) 
(B.14) 

(B.15) 



The energy and momenta of the bosons arc fixed for given values of Si and S2. The 
momentum of the W~ boson in E can be written as: 



Pw- = (Qi, Bi sin 9, 0, Bi cos 9), 



(B.16) 



where Qi and Bi are real and positive functions of s', Si, 82, cos 9 and v. The velocity 
of the l^~-boson in the laboratory system is then: 



Pi 



El 



i = 1,2. 



With 
and 



Pe- + Pe+ = Pw- + PW+ 



Pw- = Si, 
two sets of solutions may be obtained: 



Pw+ — -52, 



-Bl,2 = 
Ql,2 = 



{s' — S2 + si)-\/l — v'^ {v cos 9 ±b) 
2V^(l-v^cos'^9) ' 
{s' -S2 + si)VT^^ (1 ± bv cos 9) 



2v^(l --u^cos^e) 



(B.17) 

(B.18) 
(B.19) 

(B.20) 
(B.21) 



Here, we used the abbreviation 



l,si,s'(l - r2c()s2/?) 

y {S' -S2 + Sif{l-V^y 



1 - 



(B.22) 



The number of solutions depends on cos 6', v and b. By definition, B is real and 
positive. There is no solution, when t;cos^ < —b, one solution for |t;cos^| < b, and 
two solutions exist for v cos 9 > b. 
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With the given solutions for B and Q and eq. (B.ll), the relation between the W 
production angles in the two Lorentz systems is found: 



cos^' 



B cos 9 — vQ 



^(1 - v^)B^ sin^ e+{Bcose- vQf ' 



(B.23) 



For the limiting case of on-shell W pair production the transformation (B.23) is in 
accordance with a similar transformation given in [54]. 



C Neutral Current Kinematical Functions 

In section 2.2.3, we use two kinematical functions known from the the study of neutral 
current process [55]: 
The function gg^ is: 



(cos6l;si;s2,s) 



-(1 + COS^ 6l)^422(si;S2,s) 

1-3 cos^ 9 , X 3 , ^ ^ - , , 
H ^ si[s + S2)- 1 - 2/:(si; S2, s) 



(C.l) 



SS2 



A 4 V ' ' si — s-2 — s 

where ^422 is also known from different contexts ([9],[56]-[61]):^ 

Qa22[S] Si, S2) = C[S] Si, S2) - 2. 



S - Si - S2 



(C.2) 



The function (cos^, s, si, S2) is: 



^23?(C0S6I,S,S1,S2) = g(l + COS^6')^233(s;Sl,S2) 



3 3 



1-3 cos^ 9)s [£(si; S2, s)2s2(si - S2) + (s - si - 3s2)] 



X [C{S2] s, Si)2si(s2 - Si) + {s-S2- 3si)] , 



(C.3) 



with [37] 

^233(5; Si,S2) 



A2 
4s 



{C{s2;s,si)C{si;s2,s) 



SSi{s — Si)^ + SS2{s — S2Y + SiS2{si — S2Y 
+ (s + Si + S2) [>C(S2; S, Si)2s (s - S2)^ + Si(s + S2 

+ C{si\ S2, s)2s (s - Si)^ + S2(S + Si - 2S2 
+5s^ - 4s(si + S2) - (si - S2)^] } . 



(C.4) 



^We explicitely agree with [9, 56, 57]. 
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